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ABSTRACT 

We propose a 2-temperature radial dynamical model of plasma flow near Sgr A* and fit the 
bremsstrahlung emission to extensive quiescent X-Ray Chandra data. The model extends from sev- 
eral arcseconds to black hole (BH) gravitational radius, describing the outer accretion flow together 
with the infalling region. The model incorporates electron heat conduction, relativistic heat capacity 
£N| . of particles and feeding by stellar winds. Stellar winds from each star are considered separately as 

sources of mass, momentum and energy. Self-consistent search for the stagnation and sonic points is 
' performed. Most of gas is found to outflow from the region. The accretion rate is limited to below 1% 

of Bondi rate due to the effect of thermal conduction enhanced by entropy production in a turbulent 
flow. The X-Ray brightness profile proves too steep near the BH, thus a synchrotron self-Compton 
^sO ' point source is inferred with luminosity L ~ 3 • 10 32 erg/s. We fit the sub- mm emission from the inner 

flow, thus aiming at a single model of Sgr A* accretion suitable at any radius. 
Subject headings: accretion, accretion disks, Galaxy: center 
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1. INTRODUCTION 



Our Galaxy is k nown to host a supermassive black hole (BH) with mass M k 4.5 • 1O 6 M at a distance R k 8.4 kpc 



X 
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7^1 (|Ghez et al.l F2008). The BH exhibits very low luminosity state probably due to inefficient feeding and cooling. The BH 
Q ■ is fed by s tellar winds within s everal arcseconds from the compact object roughly around the radius of BH gravitational 
' influence (jCuadra et al.ll2008l ). The stellar winds are expelled at large spee ds. They collide, heat up to ~ 10 7 K and 
c/2 ! emit bremsstrahlung X-Rays, observed by Chandra (jBaganoff et al.l [2003h . A small fraction of mass accretes onto 
the black hole is thus producing the emission in sub-mm and othe r wavebands. H owever, the inferred accretion rate 
withi n several Sch warzschild radii is 2 orders of magnitude lower (|Marroneil2007f ) than the inferred Bondi accretion 
rate (|Bondilll952l ) at several arcseconds. This disparity is resolved in a present work with a point source revealed 
■ coincident with Sgr A*. A brief account of observations is made in §0 The dynamical model is outlined in §[31 The 
I/"} ' results are discussed in § |4j 
ITl '. 

^ , 2. OBSERVATIONS 

We analyze ~ 1 Ms of Chandra exposure of Sgr A* and central arc seconds (|Muno et al.1 12008) significantly improving 
^v^j , over the previously analyzed 41 ks exposure (Ba ganoff et al.ll2003h . As we are interested in the quiescent emission, 
we bin the observations in 628 s bins and exclude the flaring states with counts triple the mean. We also subtract 
' the point sources and extended bright emission zones, thus extracting the quiescent count surface brightness profile 
within 5". Having the extensive data we are able to perform the subpixel spatial binning in rings 0.125" thickness 
^ ■ owing to dithering of spacecraft. The counts from four 90deg ring segments centered at Sgr A* are compared in order 
to test the viability of the radial model. It appears that within 2" the counts do not differ significantly between ring 
segments, but the variation was found at > 2". The point spread function (PSF) is extracted by observing the nearby 
binary J174540.9-290014. 



3. STELLAR WINDS AND DYNAMICAL MODEL 

Feeding of the black hole should be a starting point of any accretion model. This approach helps to eliminate a 
number of arbitrary boundary conditions. A set of ~ 30 wi n d em itters is believed to supply almost all the matter 
into the feeding region of Sgr A*. Following ICuadra et al.l (|2008| ). we identify the import a nt wind emitte r s, fin d 
the wind speed s and ejection rates. We obtain the orbital data from iPaumard et "ail (|2006f ); iMartins et al.l (2007); 
ILu et all (2009), assuming the stars either belong to the disk or taken to have the minimum eccentricities. As we 
are constructing the radial model, the radial feeding function q(r) is produced by smoothing wind inputs over radius 
between the apocenter and the pericenter for each star (see Fig. QJ. The averaged wind velocity is found as a root- 
mean-square average over stars weighed with the ejection rate. We do not account for orbital velocities of stars in 
energy input as feeding is dominated by onl y a few stars close to the BH. S2 star is included into the calculation as it 
may eject more matter ([Martins et al.l f2008). than falls onto the BH. 

The dynamical model has sources of mass, radial momentum and energy due to winds ([Lamers k Cassinellil H999). 
The main feature of the model is the electron thermal conduction proportional to the temperature gradient with 
conductivity K = 0.1^//csT , e /m e n e r ([Johnson fc Quataertj|200"7t ). Here we have included the factor of 5 inhibition of 
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Fig. 1. — Mass input into the central region around the BH on the left panel. Root-mean-square wind velocity v w on the right panel. 




conductivity in tu rbulent magnetic field (|Naravan fe M ed vcdcvl 1200 li ). So computed heat flux appears to be below 
the saturated flux (|Cowie fe McK ce 1976). In reality the inner accretion flow is collisionless and the outer accretion is 
weakly collisional, however, the single prescription for conductivity gives the reasonable approximation. The realistic 
behavior of electrons in the inner flow is achieved by using their relativistic heat capacity in the equations of motion, 
which naturally leads to the ratios T p /T e up to 10 even in adiabatic flows. For completeness of the theory we add the 
direct energy transfer into electrons and protons equiva lent to the entropy production, which happens due to viscosity 
in the rotating flow or due to dissipation of turbulence (|Shcherb akov 200cf) . We assume that the fractions f e and /j of 
available gravitational energy goes to electrons and protons. The Coulomb collisions are included for numerical stability 
to balance the electron and proton temperature in the outer flow, though they do not have the significant dynamical 
effect. The effect of gas composition is accounted for by introducing the effective mass m av « 1.245m// per electron 
and correspondent reduction in ion gas p ressure. These numb ers are taken for the solar abundance of elements, which 
seems to be reasonable for stellar winds dNaiarro et al.ll2b04f ). Paczhynski-Wiita gravitational potential is employed. 

The proposed system of equations has no artificial boundary conditions, but it appears to have an unmatched 
complexity. We self-consiste ntly solve for t he positions of the stagnation point, where gas velocity is zero, and the 
inner isothermal sonic point ()Quataertll2004h . The heat flux is set to zero at the point, where dT e /dr = in the Bondi 
solution near the BH. The outer boundary is either taken to be the isothermal sonic point in the outflow or the point 
with slightly higher density (for numerical stability). The relaxation technique is used for the 2-temperature system 
between the inner boundary and the stagnation point, whereas only shooting works outside the stagnation point. 

4. RESULTS 

Having produced a bunch of dynamical models, we convert the temperature and density radia l profiles into the 
surface brightness profile. We take the up-to-date bremsstrahlung emissivities (see ([Gouldl 119801 ) and errata) and 
account for emission by heavy elements, excluding iron. We calculate the spectrum along each ray and convert it to 
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Fig. 2. — Plots of surface brightness S in observed counts per pixel squared (1 pixcl=0.5"). Blue curve with error bars shows the 
observations. Extended emission model is on the left panel: the green curve shows brightness smoothed by PSF. The model of a point 
source with L = 3 • 10 32 erg/s (red) and the residual (green) are on the right panel. 



counts by applying the solar metallicity interstellar absorption (Morrison & McCammon 1983) with hydrogen column 
Nh = 10 23 cm -2 and convolving with the response of Chandra. Then we apply the PSF blurring and compare with 
the observed surface brightness profile. The outer surface brightness profile can be reasonably fitted by a model with 
f e = 0.15, fi = 0.05, which also leads to the ratio T p /T e « 15 near the BH. However, we find that the inner part 
of surface brightness curve is too steep for any extended emission. There should be a point source in the center 
accounting for ^ to | of central surface brightness and having the unabsorbed luminosity about L = 3 • 10 32 erg/s 
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of monoenergetic 4 keV photons. About the same luminosity is expected from the synchrotron self-Compton (SSC) 
process near the BH. The search for the best model with the point source co ntinues. We want the best model to 
reproduce t he observed Faraday rot ation measure RM « 50cm- 2 |M arron i2007t) and optically thick flux F R = 1.73 Jy 
at 86 GHz ([Krichbaum et all |2006[ ) . An order of magnitude consistency is achieved on the way. The accretion rate 
appears to be self-consistently limited to < 1% of Bondi value, thus the connection between the inner accretion flow 
and the outer accretion flow is established. The future versions will include the angular momentum and use X-Ray 
spectral information. 

The author is grateful to Mikhail Medvedev and Ramesh Narayan for fruitful discussions, Daniel Wang and Feng 
Yuan for comments. 
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